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Abstract 


In an earlier article [3], we presented an algorithm that can be used to rigorously 
check whether a specific cosine or sine polynomial is nonnegative in a given interval 
or not. The algorithm proves to be an indispensable tool in establishing some recent 
results on nonnegative trigonometric polynomials. See, for example, a and [5]. It 
continues to play an essential role in several ongoing projects. 

The algorithm, however, cannot handle general trigonometric polynomials that in¬ 
volve both cosine and sine terms. Some ad hoc methods to deal with such polynomials 
have been suggested in [3], but none are, in general, satisfactory. 

This note supplements [3] by presenting an algorithm applicable to all general 
trigonometric polynomials. It is based on the classical Sturm Theorem, just like the 
earlier algorithm. 

A couple of the references in [5] are also updated. 
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1 A Summary of the Algorithm Proposed in |[3| 

Cosine Polynomials 


1. Any cosine polynomial C{x) = cos{kx) can be rewritten as an algebraic poly¬ 

nomial P{X) in the variable X = cos(x). The nonnegativity of C{x) in a given 
interval follows from that of P{X) with A in a corresponding interval. 

2. The nonnegativity of P{X) can be rigorously verified using the classical Sturm The¬ 
orem (see m)- 


With the help of the symbolic manipulation software MAPLE, these steps can be auto¬ 
mated with simple-to-use commands. The same is true for the algorithm involving sine 
polynomials. 


Sine Polynomials 

1. Any sine polynomial S{x) = X]fc=i rewritten as a product of sin(x) 

and an algebraic polynomial P{X) in X. The nonnegativity of S{x) can be deduced 
from studying the signs of P{X) in suitable subintervals of a corresponding interval. 

2. The signs of P{X) can be determined rigorously using the classical Sturm Theorem. 


General Trigonometric Polynomials 


Some ad hoc methods are proposed, but none are satisfactory. 


2 New Algorithm for General Trigonometric Polynomials 


The new procedure is based on the elementary identities 


2T 

sm(x) = 


cos(x) 


1 -T2 

TTt2’ 


T = tan (I) . 


( 1 ) 


1. Any trigonometric polynomial can be rewritten as Pi (A)-|-sin(x) P 2 (A), where Pi (A) 
and P 2 (A) are algebraic polynomials in A. 

2. Applying ([1]) results in a rational expression, the numerator of which is an algebraic 
polynomial in T and the denominator is a power of (1 -|- T^). 

The nonnegativity of the original trigonometric polynomial in a given interval follows 
from that of the numerator of the rational expression in a corresponding interval. 
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3. The nonnegativity of the numerator can be rigorously verified using the classical 
Sturm Theorem. 


Remark. Although the new algorithm can also be used for a purely cosine polynomial or 
a purely sine polynomial, it produces an algebraic polynomial with degree twice as that 
produced by the older algorithm. Hence, for purely cosine or sine polynomials, the older 
algorithm is preferred. 


3 Examples and MAPLE Commands 


Example 1. Let us illustrate the new procedure by proving that 

. / N 3 , , , , sin(2x) cos(2x) . , 

A(x) = - + sm(x) + cos(x) H-- —--\ --—^ > 0, a:€[0, vr]. (1) 

0 A A 

A{x) can be rewritten as 

— + cos (x) + sin (x) + cos^ (x) + sin (x) cos (x). 

After applying ([T]), it becomes 

- IST^+ 40r + 21 
10(l + r2)2 

As X varies from 0 to tt, T varies from tan(O) to tan(7r/2), or from 0 to oo. Thus, ([T]) is 
true if we can show that the numerator of ([2]) is nonnegative for T E [0, oo). The latter fact 
can be established by showing that the numerator, which is a polynomial, has no root in 
[0,oo), using the classical Sturm Theorem. 

□ 


Using MAPLE 


For more complicated trigonometric polynomials, the calculations required in the two steps 
are often very time-consuming and error-prone. They can be automated using MAPLE. 

The following command defines a MAPLE function that takes a trigonometric polyno¬ 
mial as input and outputs the desired algebraic polynomial in T. 

pt := numer(subs(sin(x) = 2*T/(1+T~2),cos(x) = (1-T“2)/(1+T“2), 
expand(f))); 
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After issuing the above command, the next two affirm the assertion of Example 1. 


p := pt(3/5 + sin(x) + cos(x) + sin(2*x)/2 + cos(2*x)/2); 

Sturm (p, T, 0, infinity); 

The first command produces the numerator of the fraction in ([2]) and assign it to the 
variable p. The second command uses the MAPLE built-in command sturm to determine 
the number of roots of pi in the specified interval T G [0, oo). In this example, the output 
from the second command is 0. 


sinf2x^ cosf2x^ 

Example 2. Find the minimum of sin(a:) -|- cos(x) H- - -1- - -, x G [0, tt]. 

The problem is equivalent to the determination of the least value a such that 

„ , , . / s / N sin(2x) cos(2x) 

Ba{x) = a + sm(x) -|- cos(x) H -^ —- H- ^—- >0 x G [0, tt]. 

Applying pt to Ba{x) produces the polynomial 

P 2 {T) = {2a - l)r^ -h (4a - 6)T^ + 8T + (2a + 3). 

Thus, we need to solve the equivalent problem of determining the least a such that P 2 {T) > 0 
in [0,oo). It is well-known that a necessary condition for such an a is to be a root of the 
discriminant of p 2 {T). Using MAPLE, the discriminant is found to be 

A = 16384(2a - l)(2a - 3)(8a2 -F 12a - 8). 

It is then easy to verify that one of the four roots of A, namely, a = 3(-\/3 — l)/4, is the 
desired answer. 

□ 

Example 3. Find the minimum of sin(x)-|-cos(x)-t-sin(2x)-|-cos(2x), x G [0, vr]. The same 
arguments as in the previous example prove that the minimum is a root of the equation 
(the lefthand side is the discriminant) 

32768 - 8192 - 211968 -k 165888 a + 27648 = 0. (3) 


However, the solutions no longer have simple representations. Numerical computation gives 
the minimum as 

a « 1.040168473. 


□ 
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Example 4. In a recent project, we prove that 


sin(x) sin(2x) , , 23 , 

2!^ + + 2^ + sm{4,.) + ^ . 5, 

(which are the first two cases of the more general polynomial) 

sm(kx) 23 . , 

-TTT + 7^ • (^ +1)’ ?^ = 3,4, •••, 


k=l 


n-k + 1 125 


are all nonnegative in [0,oo). Let us apply our new procedure to the first two polynomials. 
The corresponding algebraic polynomials are 

276 T® + 1750 + 828 - 7000 + 828 + 3250 T + 276 

and 

138 T® - 875 + 552 T® + 7375 T® + 828 - 9025 T® + 552 + 1925 T + 138 

respectively. In both cases, the Sturm Theorem shows that the polynomials are nonnegative 
for T E [0, oo). 

□ 
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